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Abstract 

We discuss certain specific features of the calculation of the critical temperature of a dilute 
repulsive Bose gas. Interactions modify the critical temperature in two different ways. First, 
for gases in traps, temperature shifts are introduced by a change of the density profile, arising 
itself from a modification of the equation of state of the gas (reduced compressibility); these 
shifts can be calculated simply within mean field theory. Second, even in the absence of a 
trapping potential (homogeneous gas in a box), temperature shifts are introduced by the 
interactions; they arise from the correlations introduced in the gas, and thus lie inherently 
beyond mean field theory - in fact, their evaluation requires more elaborate, non-perturbative, 
calculations. One illustration of this non-perturbative character is provided by the solution 
of self-consistent equations, which relate together non-linearly the various energy shifts of the 
single particle levels k. These equations predict that repulsive interactions shift the critical 
temperature (at constant density) by an amount which is positive, and simply proportional 
to the scattering length a; nevertheless, the numerical coefficient is difficult to compute. 
Physically, the increase of the temperature can be interpreted in terms of the reduced density 
fluctuations introduced by the repulsive interactions, which facilitate the propagation of large 
exchange cycles across the sample. 



Introduction 



The calculation of the effects of interactions in dilute gases is often considered as a classical 
textbook problem, which is generally solved with the help of cluster techniques in statistical 
mechanics; well known examples are for instance the second virial correction to the pressure, 
the heat capacity, or to the magnetic susceptibility for a gas of spin 1/2 particles E]- For 
gases at low temperatures, all these corrections are expressed in terms of a single parameter, 
the scattering length a, so that the result is simple. At first sight, the calculation of the first 
correction to the critical Bose-Einstein condensation (BEC) temperature in a dilute repulsive 
Bose gas seems to raise a similar problem. Nevertheless, its solution was not well understood until 
recently, as illustrated by a large collection of contradictory results in the literature (see |3] and 
references contained). Several reasons underlie the confusion. First, for a homogeneous gas at low 
temperatures, it turns out that a mean field treatment of the interactions leads to an exactly 
zero shift of the critical temperature (the critical value of the chemical potential is shifted, but 
this does not affect the critical temperature); combining mean field theory with additional - and 
not necessarily well controlled - approximations can then lead to results which are approximation 
dependent, providing arbitrary a power variations and even sign. The second reason is that the 
problem is actually not as simple as it looks. It is essentially non-perturbative at several levels 
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3 4, even if the leading correction term which emerges from the calculations ends up being simply 
linear in a; we discuss below in more detail this peculiar feature of the calculations. In retrospect, 
this difficulty is actually not so surprising, since the properties of a Bose gas are non-analytic 
around a = 0, where the system is at the border of stability (a Bose gas collapses at condensation 
as soon as a becomes negative). The linearity in a is therefore far from obvious, and indeed it 
turns out that the next correction |2j is not simply ~ a 2 but ~ a 2 log a, clearly manifesting the 
non-analytic character of the problem. 

The purpose of the present article is to return to the considerations of [j^ and 0j, with more 
detailed discussion of some of their physical aspects. In particular, we wish to emphasize that the 
linearity in a of the leading term in the correction to the critical temperature may be missed if 
a non self-consistent theory is used to calculate the effect; indeed, some calculations found in the 
recent literature and based on non self-consistent models provide different results, first corrections 
~ a log a for instance. But, before we come to this point, in order to avoid any confusion, we 
carefully distinguish between two effects which are discussed in this context; the former takes 
place in traps only, the latter also in a uniform gas contained in a box. 

1 Compressibility effects in a trap (mean field) 

For a Bose gas of N atoms contained in a trap at a given temperature, Bose-Einstein condensation 
first sets in at the point where the density is maximum. Repulsive interactions tend to make 
the density more uniform, and hence will in general lower the density at this point, so that 
they will decrease the transition temperature with respect to an ideal gas. Several authors have 
studied this effect; for a review, see ref. 6 (§ V-B) and references therein. Here, for the sake 
of simplicity, we limit ourselves to the discussion of isotropic traps and to the thermodynamic 
limit 1 . This excludes finite size effects, which are also discussed in 6 ; a finite size system does 
not undergo a sharp phase transition and thus the critical temperature, and its shift, depend 
on a somewhat arbitrary definition of the critical temperature 2 . In the thermodynamic limit, 
this arbitrary character disappears; the width of the trap becomes large so that the effects of 
the trapping potential may be treated within a semi-classical approximation. Around the center 
of the trap, one can then completely ignore the effects of the potential, and the condensation is 
obtained when, locally, the critical conditions are reached. This happens when the local degeneracy 
parameter, the product n(0)A 3 , is equal to the critical value of this parameter in a homogeneous 
system, where n(0) is the number density at the center of the trap and A = h/ \f2-KmkBT the 
thermal wavelength. For an ideal gas, this value is (3/2 — 2.61..; for an non ideal gas, it is slightly 
shifted, see §|21 The question now is to relate ro(0) to N in order to obtain its critical value as a 
function of the temperature, or conversely the critical temperature as a function of N. 

1.1 Equation of state 

In a gas in equilibrium, the local density of the gas adjusts locally so that the pressure gradient in 
the gas compensates exactly for the force exerted by the external potential. The equation of state 
of the gas therefore plays an essential role in the determination of the density profile in a trap. A 

1 The thermodynamic limit in a trap is denned as the limit where the number N of particles tends to infinity, 
the frequency ui of the trap goes to zero, and the product Nuj 3 remains constant. 

2 In an ideal gas, the number of particles in excited states N e can not exceed a maximum value ./V™ ax , obtained 
when the chemical potential is equal to the single particle ground state energy. The usual definition of the critical 
value N c for the total number of particles TV in a finite trap is N c = 7V™ ax . This convention is convenient, but 
remains somewhat arbitrary, since it relates the value of one physical quantity, N, to the upper limit of another 
physical quantity N e . It therefore conveys the idea of N e increasing until it reaches its absolute upper limit, and then 
saturating. In reality, a large number of particles has already accumulated in the ground state when N = 7V™ ax , 
and N e actually never reaches this upper limit (except in the limit N — > 00). 

Another definition of the critical point may be obtained by plotting the number No of particles in the lowest 
state as a function of temperature (at constant N), and taking the inflexion point to define the critical temperature. 
This is in a sense more physical, but still arbitrary, and this different convention leads to different values of the 
finite size effects. 
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Figure 1: Figure (a) shows the number density n in the gas as a function of the chemical potential 
/i, at constant temperature T (n, measured in unities of A -3 , is dimensionless). This curve is 
obtained within simple mean field theory; the critical value of /x is /z c ; beyond this value, the slope 
of the curve is inversely proportional to the scattering length a. Figure (b) shows the pressure P 
in the gas, which is the integral of n over fi . Figure (c) shows the equation of state relating P to 
n, after elimination of /i; the critical value of the density is n c \ at this point, the compressibility 
of the gas is proportional to 1/a. Dashed lines refer to the ideal gas (a = 0). 



good approximation is provided by mean field theory, where the local density of the gas is written 
as: 

n = ^3 33/2 = exp^O- Afi)] +n (1) 

with: 

0Afi = Aa\ 2 n (2) 

where f3 = 1/ksT is the inverse temperature, 33/2(2) the usual Bose function |2] and no the 
density of condensed particles. Figure 1-a shows the variations of the density 3 as a function of the 
chemical potential fi, assuming that the gas is repulsive (a > 0) - see for instance figure 10 of ref. 
[7] for a geometrical construction of this curve. Integration over the chemical potential provides 
the pressure (fig. 1-b); eliminating fi leads to the equation of state in mean field approximation 
(fig. 1-c). For comparison, the curves are also shown for the ideal gas (a = 0). As expected, 
repulsive interactions tend to increase the pressure of the gas and to reduce its compressibility. 



1.2 Total number of particles in a trap vs. density at the center 



In a trap that is sufficiently large (thermodynamic limit), the effects of the external potential can 
be treated semi-classically, which allows one to extend equations (11121) to inhomogeneous systems. 
We now assume that the gas is not Bose condensed and define the effective chemical potential /x e ff 
by: 

fi cS (r) = /x - V(r) - 4a\ 2 n(r)/(3 (3) 

and find: 



n(r)X 3 - 33/2 [z = e^ ff(r) 
But equation (0} can be inverted as: 

/icff = /3" 1 J(nA 3 ) 

where J is the logarithm of the inverse function of 33/2- By differentiation, one obtains: 

VMeff = /3 _1 A 3 J'(?iA 3 ) v« 

3 We take this opportunity to note that fig. 2 of |3| is inaccurate: at the critical value fj, c of the chemical potential, 
the curve should not have a kink, but a continuous slope. 



(4) 



(5) 



(6) 
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where J is the derivative of J. Combined with (J3J, this equation yields: 



-13 V V(r) = [A 3 J'{n\ 3 ) + 4aA 2 ] V «(r) (7) 

which relates the local variations of the potential and of the density. In general, the term in 4aA 2 
on the right side of this equation is a small correction to the term in J', so that the relation 
between the two gradients depends only weakly on a. Nevertheless, in regions of space where the 
gas is almost condensed, the derivative J becomes very small and, for a given \/V, the density 
gradient \/n becomes much larger and strongly a dependent; this phenomenon can take place at 
the center of the trap, as we discuss in more detail in § 

Equation (7J allows one to calculate n(r) at every point of the gas as a function of n(0) and, 
through an r integration, of the total number of particles A. At the critical point, n(0) is fixed by 
the condition n(0)A 3 = C3/2 m the mean field approximation 4 , so that this calculation provides a 
relation between the critical values of A and T. For a harmonic trap, where the size of the single 
particle ground state is af lo = y/h/rnuj, Giorgini et al. [H] show that, to lowest order, this critical 
temperature is shifted with respect to its ideal gas value T® by an amount 5T C given by: 

^ = -1.3— (8) 

J- c Ciho 

Since TV 1 / 6 /ciho ~ (Aw 3 ) 1 / 6 , the relative change of temperature remains finite in the thermo- 
dynamic limit. Significant critical temperature shifts can indeed be observed experimentally; for 
instance, recent experiments by Gerbier et al. [U] have reported observations of relative shifts as 
large as 10% in a Rb gas, in good agreement with JHJ; early Monte Carlo calculations on a trapped 
Bose gas also predicted a temperature shift dominated by the mean field ^01 ■ 
A few remarks may be useful at this stage: 

(i) the shift 5T C describes the case where the total number of atoms A is kept constant; if, 
instead, the density of the gas at the center of the trap were kept constant, the shift would vanish 
(within mean field theory). Within a factor, A and n(0) are conjugate variables 5 , the former 
extensive, the latter intensive, somewhat similar to energy and temperature in thermodynamics. 
Neither of them is measured directly in experiments, since optical measurements provide the 
"column density" (density integrated along the line of sight), an intermediate physical quantity. 
Conceptually, there is nevertheless a preference for using n(0) as the relevant variable: while the 
calculation of the shift 5T C at constant A gives different results for all possible forms of traps 
(parabolic, quartic, or even more complicated as that shown schematically in figure 2), in the 
thermodynamic limit, all lead to the same value of n(0); this result is clearly more universal. 

(ii) the physics involved in ST C is only remotely related to the physics of the Bose-Einstein 
transition; it is more closely associated with the effects of interactions on the density profile in 
trap, which is in turn determined by the equation of state. In other words, the calculation of 
A involves many atoms which play actually no role in the condensation, because they are in 
regions of space where the density is too low, almost or completely in a classical regime. This is 
particularly obvious for the trap of figure 2 where, clearly, all the atoms in the external part of 
the trap are counted in A, with a weight proportional to the volume of this part, while they play 
no role whatsoever in the condensation phenomenon. 



2 Correlation effects in a uniform gas 

We now come to the central part of this article and discuss another shift of the critical temperature 
of totally different physical origin. This effect is not related to an external trapping potential but 

4 In all this article, we assume that the interactions can be described by a single parameter a, and therefore have 
no k dependence. In this case, one can easily show |3||7| that, within mean field approximation, the critical value 
of the degeneracy parameter is unaffected by the interactions. 

5 When the potential is treated semi-classically, one can show that N = ff^^O 2 ) x 8LogZ/dn(0) where Z 
is the partition function with the usual notation z = exp f3/i. 
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Figure 2: Example of a shape of a trap where only the central part plays a role in Bose-Einstcin 
condensation; all the atoms in the other part are counted in the total number of particles N but 
do not take part in the phenomenon. 



emerges from correlations between the particles; it is actually most conveniently calculated for 
a uniform gas contained in a box. As recalled in the introduction, the study of this shift has 
a long history of contradictory results, predicting either an increase or a decrease of the critical 
temperature, as well as variations with different powers of a (see [H]). It is now understood that 
correlation effects taking place in the non condensed gas just above transition lead to an increase 
of the critical temperature T c , which is proportional to a in leading order: 

^ - c onV» (9) 

where c is a numerical positive coefficient. In linearity emerges as an exact result from a 

scaling analysis of the complete perturbation series. 

In this text, we give another version of the argument leading to this linearity, with more em- 
phasis on simple properties of non-linear self-consistent equations, and on the physical mechanism 
they contain. Behind the increase of T c is a modification of the effective energies that determine 
the populations of the particles with low momenta. Within the framework of the non-linear self- 
consistent equations introduced in and without elaborate calculations, we summarize the 
basic ideas that lead to the shift of the transition temperature, and comment on various points 
involved. 

2.1 Non-linear equation for the single particle energy 

In a non condensed gas, the diagonal elements of the single particle density matrix can be 
written as: 

' s /30k-M+£ k ) _ i\ ~ 



Pk = 



(10) 



where (3 is the inverse temperature, fj, the chemical potential and the effective energy is the sum 
of the free particle kinetic energy e k = ft 2 k 2 /2m and the energy shift introduced by the 
interactions. As in [7], we write 6 : 

/3£ k = 4a\ 2 n-8(j) 2 (J^j J d 3 k p k , J d 3 q p k+q pv_ q (11) 

This equation is not exact, but a self-consistent "one bubble approximation" to £ k ; it neverthe- 
less allows a qualitative discussion of the properties of the more general theory and provides a 

6 In Ursell theory, equations 1101 and 1111 are naturally obtained in a self consistent second order approximation. 
Note that the energy shifts are then distinct from the self energies defined in Green's function theory; in 
particular, they provide no information on the evolution of the quasiparticles, but just the static populations p k . 
See ref. 1121 for a discussion of the appearance of exponentials exp/3Sk in Ursell operator theory. 

In Green's function formalism, equations I1UI and 1111 can also be obtained with the help of additional approx- 
imations. For instance, one can neglect of the Matsubara frequency dependence of the self energies, which allows 
one to reconstruct in p k the Bose-Einstein distribution function appearing in 1101 . and recover the same equations. 
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reasonable estimate of the leading correction for AT C in the homogeneous system. The first term 
on the right side of is simply the mean field term; in the low temperature regime where the 
interactions are described in terms of the scattering length a only, it depends only on the density: 

W Pk (12) 
Because this mean field is independent of k, it has no effect whatsoever on the critical density 
(at a given temperature). The second term is more interesting; it corresponds to the effect of 
correlations due to interactions in the gas, and introduces a k dependence of the energy shift Sk- 
We will see that this shift tends to "harden" the free particle spectrum around k = 0. 
The condition for Bose-Einstein condensation reads: 

H = S (13) 

so that, in terms of the "spectrum" at the critical point: 

W(k) = p [e k + S k - S ] (14) 

eqs Hl()(l and (|ll|l provide the following self-consistent equation: 



W(y) -**-»(5)'(s)7 **'/*« 



1 1 



e W(k') _ 1 e W(k'-q) _ l 



1 1 



e W(k+q) _ 1 e W(q) _ l 



(15) 



From the spectrum we can obtain the populations and the critical density n c . Actually, what we 
wish to obtain is the change An c of the critical density (at constant temperature) with respect to 
the ideal gas value n° = C3/2A 3 ■ We assume that only small values of k contribute to this change 
(we discuss why this is true below); since W(k — 0) vanishes, we can replace [exp W(k) - If 1 by 
1/W(k) to obtain: 



An c 



d 3 k 



/3e k - W(k) 



Pe k W(k) 



(16) 



to leading order. Note that, as mentioned above, the mean-field term does not enter the problem 
anymore, and therefore does not influence the critical density or temperature. 

At this point, it is easy to convert the change of density l|16f) at constant temperature into a 
change of temperature at constant density. The critical condition in the presence of interactions 
can be written as: 

n h 

nX3 = <o ^ ^3/2 = C3/2 + A 3 A?i c (a, T) (17) 

where An c (a,T) is given by (|16fl . Equation l|17|) defines the new critical line in the temperature- 
density plane (fig. 3). Any point on the initial critical line (ideal gas) can be moved by a first 
order change of either the density, or the temperature, and reach a point on the new critical line; 
the only condition is that the change of the product nA 3 should be equal to A 3 An c (a,T) to first 
order 7 . But the same change of nX 3 can be obtained, either when An = with a change of 
temperature AT , or when AT = with a change of density An. In fact, for any motion in the 
temperature-density plane, the relative variation of the left side of 1|17|) is given by: 

A(nA 3 ) _ An 3 AT 
(nA 3 ) ~ ~ ~ 2~ ( ' 

The same relative variation can therefore be obtained, either with a relative change An/n at 
constant T, or with a relative change AT/T at constant density, provided the ratio between the 

7 To first order, we can ignore the difference between An c (a, T + AT) and An c (a, T). 
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Figure 3: Critical line in the density-temperature plane, for the ideal gas (upper line) and for the 
interacting gas (lower line). Starting from a point on the first critical line, arrows indicate how 
a change of density at constant temperature, or conversely, move the point to the new line, as 
discussed in the text. 

two changes is —2/3. Finally, the leading term of the change of critical temperature AT C is given 
by: 

AT C 2 r d 3 k 

We will see in § 12.1.21 that the change of the spectrum at the critical point manifests itself 
primarily in a small k region centered at the origin and of width Afc given by: 

Ak = k c ~-^ (20) 

In this region, the correction Ek — Eo is comparable to the free particle energy ek, or may even 
dominate in VF(k); this property is also discussed in detail in [31 EI • The physical origin of the 
change of energy is simple: the atoms with very low velocities are extremely sensitive to even 
very small effects of the interaction potential. They can therefore rearrange themselves in space to 
minimize their repulsive energy, reaching a smaller value than the mean field prediction. On the 
other hand, the atoms with higher k values than Ak have too much kinetic energy to do so. 

This rearranging process can be regarded as an analogue, in terms of interaction potential and 
correlations, of the effect of a weak external potential on particle positions discussed by Lamb et al. 
[T3*| . These authors point out that the character of the Bose-Einstein condensation is qualitatively 
modified even by a weak external potential, transferring it from momentum space to real space. 
Here we have another illustration of the extreme sensitivity of a gas just above transition to any 
potential, but in the space of relative positions instead of ordinary position of the particles. 

2.1.1 Non self-consistent approximation: spurious logarithms 

The non- linear equations 1)10111(1 . or l)15|l. are not easy to solve. One can start by assuming that 
Ek provides only a small correction to the free particle spectrum, so that the right side of (flip 
can be calculated with Ek = 0, which provides a first approximation for Ek - this operation can 
be regarded as a first iteration of an infinite process leading to a fully self-consistent solution. 
The critical condition (|13f) then gives an implicit equation in /i to determine the critical value of 
the chemical potential (still within this first iteration). For this value of /x, one can then use (f 1 Of) 
to calculate the new populations /5k, and finally n by integration over k. This is precisely what 
was done in for values of a/ A ranging from 10 -3 to 10 -2 , leading to the following value of the 
coefficient c (obtained by dividing the relative change of T c by an 1 / 3 ): 

c~0.7 (21) 

One can of course also make two (or more) iterations: again choose /i, calculate the first 
approximation for Ek, the first approximation for pk, inject them into (|11J) to obtain the second 
approximation for Ek, and then only determine fj, by the condensation condition. The second 



/?e k - W(k) 



f3e k W(k) 



(19) 
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approximation for eventually provides the critical density by integration over k. This procedure 
leads to a progressive hardening of the spectrum, as illustrated by Fig. 15 of [7], and typically to 
higher values of c than suggested by the first iteration approximation, eq. (|21|) . 

But what is the validity of these calculations? Knowing that the presently accepted value is 
c ~ 1.3 [HUH, at first sight 

one could consider this calculation a success, especially in view of 
its simplicity. However, success may be pure coincidence! In fact, an analytical study of the single 
iteration procedure shows that it does not lead to a linear dependence of AT C in a, but actually 
contains spurious logarithms in a; these logarithms will dominate in the limit a — > and destroy 
the linearity - we come back to this point in § 12.2.11 but see also § 5.1 of In other words, if 
the numerical calculations leading to 1(211) had been made with a different range for the parameter 
a I A, a different value for the coefficient c would have been obtained, as illustrated by the table of 
§ 5.1 in The validity of the first iteration leading to eq. I(21|) is therefore difficult to assess a 
priori; see |16j for another example of a calculation which predicts a non-linear leading correction, 
an a log a in this particular case. 



2.1.2 Self-consistent treatment: linearity 

We now come back to the full equation (|15fl and discuss the properties of the self-consistent solution 
W (k) . We define the dimcnsionless function v(x) of the dimensionless variable x — kX 2 /a by: 



Inserting this into 1(15(1 provides: 



V( - X) _ 4tt (2^) 6 



W(k) 



dV 



d 3 y 



kX 2 



(a/A) 2 



(22) 



VA) 2 



(a/A) 5 



e (a/A)Mx') - 1 e ( a / A )M* 
(a/A) 2 



e (o/A)Mx+y) - 1 e WA)My) _ 1 



(23) 



If we take the limit a/A — > of the integrand on the right side of this equation (an operation 
justified below), we obtain the parameter-free self-consistent equation: 



V{X) = 



(2tt) £ 



dV / d 3 y 



1 



u(x')u(x' 



1 



1 



v(x + y) v(y) 



(24) 



(one can easily check that no other choice for the powers of a and A in the scaling factors of l|22[l 
would allow the same complete disappearance of all parameters from the equation). 

We now assume 8 the existence of a solution to 1(24(1 that, for values of x 1, is dominated 
by the free spectrum term x 2 /47r. We note that the integral is made ultraviolet (UV) convergent 
by the difference that appears in the brackets. As for infrared (IR) convergence, for small x's the 
self-consistent solution will automatically have a harder variation than x 2 (free spectrum) in order 
to make the integral convergent, a variation in x 3 ^ 2 for instance 9 ^7j. Being convergent, the value 
of the integral can be obtained with arbitrary accuracy from a finite domain of the dimensionless 
variable x, corresponding to |x| < D (where D is a pure number). Coming back to l|23|l . we see 
that the exponent in the denominator is smaller or equal to (a/ X) 2 v(D), which tends to zero when 
a/A — > (since v(D) is a pure number, independent of a). This, in retrospect, justifies the lowest 
order expansion of the exponentials inside the integrand which led us to l|24() . at least for values 
of the current variable x which are not too large (see below) . 

Finally, inserting 1)22(1 into 1)16(1 provides: 



2 

7TA 2 



dk 



(a/X) 2 v{k\ 2 /a) - (fcA) 2 /47r 
(a/X) 2 v(kX 2 /a) 



(25) 



8 This assumption is supported by numerical and analytical calculations of 0. 

9 A simple power counting argument applied to equation 1241 predicts a x 3 / 2 dependence of v(x) for small values 
of x. 
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or, with a variable change x = k\ 2 /a: 




dx 



v(x) ~ x 2 /Ait 



(26) 



v{x) 



This result is equivalent to l|16[l and shows that the density shift at constant temperature is indeed 
linear in a; consequently, the same property is true for the temperature shift at constant density. 

The values of the dimensionless variable x that contribute to the integral are smaller or com- 
parable to D. In terms of the initial momentum variables appearing in l|15|l . this domain has a 
range given by some number multiplying k c defined in (J2DJ. Going from (|23[) to 1|24|) actually also 
requires that the current variable x not be too large, namely that v(x) <ti (A/a) 2 . Since, for large 
x, v(x) is dominated by the kinetic energy term ~ x 2 , this condition corresponds to x X/a and 
therefore to kX <C 1. In other words, the values of £k obtained from l|24() are correct, but only for 
momenta significantly smaller that the maximum thermal momentum. This caveat is irrelevant 
for the calculation of the leading term of the density shift. As we will see in more detail in § 12.3.11 
what is important is the crossover value of k at which the function v(x) switches from a behavior 
dominated by the integral to a behavior dominated by the quadratic kinetic energy. From (|24|) . 
we know that this phenomenon takes place for some finite value of x which corresponds, in terms 
of k, to a multiple of k c . Finally, it is easy to see in 12U|) that, for sufficiently small values of a, the 
product k c X remains much smaller than 1, which means that only particles with low momentum 
have a role in the change of critical density. 

The conclusion is that the non-linear self-consistent equation for the effective energies leads 
to linear a dependences of the temperature shift, but only if it is treated self-consistently; other- 
wise, more complicated spurious variations are introduced - see also ref. |18| for other numerical 
illustrations of this property. 

2.2 Generalization 

We now go beyond the simple approximation (|llfl for the effective energy and generalize the 
considerations of the preceding section. In the complete theory [3] 0], the right side of this 
equation contains an infinite sum of integrals, corresponding to all diagrams in the theory. They 
start from the second order term in a/X already contained in and continue with terms of 

all orders in a/X, containing higher dimensional integrals with more elaborate structure. The 
problem then becomes significantly more complicated; we will see, nevertheless, that most of the 
conclusions of the preceding section remain valid. We first discuss the general relation between 
infrared convergence of the theory and the linearity in a of the leading term in the correction, then 
show that the linearity of the simple model remains valid within the general theory, and finally 
emphasize the difficulties arising in the precise calculation of the linear coefficient c. 

2.2.1 Linearity and infrared convergence 

It is well known that the perturbative treatment of second order phase transitions gives rise to 
infrared (IR) divergences. Here, we discuss the close relation between these divergences and a 
possible breakdown of the linearity in a of the critical temperature shift. 

Let us first come back to the single iteration approximation discussed in § 12.1.11 If we introduce 
again the effective chemical potential fj, e s = fx — 4aA 2 n//3, we see that this approximation involves 
integrals - right side of (| 1 1 11 - which, because they contain the ideal gas spectrum ek ~ k 2 , are IR 
divergent if fi c s — 0. More precisely, when /z e ff — > 0, one obtains /3E — 4aA 2 n ~ (a/X) 2 log(— /3/i e ff), 
so that the critical condition (|13fl leads to logarithms in the value of the critical chemical potential, 
— /3/ieff ~ ( a /^) 2 l°g(-V a )- But, for small k, the integrand can be written as: 



[/3(e k - /ieff)]" 1 = [(h 2 k 2 /2m) - /i eff ] 1 ~ [(fc() 2 + 1 



l 



(27) 
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Figure 4: Graphical construction to obtain the mean field correlation length, £; the lowest full line 
corresponds to the absence of IR divergences, the upper full line to an IR diverging function. The 
slope of the dashed line is proportional to (A/a) 2 and diverges in the limit a — > 0. 



where the mean field correlation length £ is defined by — fi e s — h 2 /2m( 2 , as in refs. We 
then see that ( plays the role of a scaling factor for the k dependence of the energies: 

W(k) ~ ^v(kC), (28) 

where v(x) is a dimensionless function. We can finally use this result to calculate the critical 
density shift, eq. (fTfi)l . and obtain An c oc which is proportional to (a/X 2 )y/\og(X/a). 

For extremely small values of a, the logarithmic dependence on a of the correlation length then 
dominates the shift of the critical temperature. This explains why the table of |3] gives values 
of c that depend slightly on the ratio a/A. Nevertheless, this logarithmic scaling is spurious, an 
artifact of the truncation of the self-consistent equations at the first iteration; as we have seen in 
§ 12.1.21 the complete self-consistent solution leads to a shift of T c linear in a. 

More generally, the relation between IR convergence and linearity is easy to understand. The 
critical condition (|13fl is the key to the calculation of the scaling factor £ c at the critical point, 
which in turn determines the critical density. Since Sq can be written as: 



E = — 1(0 (29) 

where /(£) is a positive integral depending on the chemical potential through £, the critical con- 
dition (fH!)l reads: 

-«=44h?) 2/<<) < 3 °> 

Graphically, as shown in figure 4, the critical value of £~ 2 is obtained by intersecting a straight 
line with large slope (A/a) 2 with the function / , considered as a function of C~ 2 - The intersection 
point is close to the vertical axis. If / has a finite limit Iq when £ -2 — > 0, we immediately obtain 
£~ 2 ~ (3Io a 2 /A 4 , and from this simple scaling factor linearity in a follows. On the other hand, if 
this function diverges when £~ 2 — > 0, the scaling factor has a more complicated a dependence. 

Ref. discusses various models illustrating the relation between the absence of IR divergences 
and the linearity in a of the leading term in the correction. For instance, the "bubble sum" , or the 
"ladder sum", both contain integrands with fractions that automatically control the divergences, 
and therefore lead to linear dependences of the density shift in a. This property explains why ref. 
|19j predicts a linear a dependence, in an approach based on the many body T matrix approxima- 
tion, similar to a ladder sum approximation. Clearly, these results do not establish the linearity in 
a for the full theory, since there is no special reason why only bubbles or ladder should be retained 
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from the full series. The disappearance of the IR divergences from the bubble or ladder sums is a 
special property of these subseries. For the reasons discussed above, one could expect the density 
corrections to be dominated by other diagrams containing IR divergences, therefore providing a 
result in a In a for instance (this is indeed the result obtained later in ref. and making the 

contribution of isolated ladders negligible. We now discuss the full theory to examine why the 
linearity in a is actually preserved. 

2.2.2 Linearity in the full theory 

The reasoning of § !2.1.2l can be extended to the more complicated structure of the effective energies 
in the complete theory. We refer the reader to the discussion given in § 4 of ref. [3]; here, we give 
only a brief discussion of the argument leading to linearity. The various terms in Ek, corresponding 
to a series of diagrams with increasing powers of a/A, may be obtained by recurrence from one order 
to the next. At each step, one includes one more power of a/A, one more k integral with a factor 
A 3 , and two more Pk's, introducing two fractions with exponentials of W(k) in the denominator. 
Now, if we make the same change of variable x — k\ 2 /a as before, we see that the change of 
integration variable introduces an additional factor (a/A 2 ) 3 , which comes in addition to an a/A 2 
from the two coefficients above; altogether, we get precisely the factor (a/A 2 ) 4 that is necessary 
to add a factor (a/A 2 ) 2 to the numerator of each fraction, exactly as in 123|) . As a consequence, 
we obtain again a parameter free non-linear equation, generalizing l|24|l . 

This new integral is more ultraviolet convergent than in the right side of 124|) since, for large 
x, two additional factors of x 2 overcompensate the additional integration factor d 3 x. As for IR 
convergence, again it is obtained because the solution v(x) has to adapt its small x behavior to 
ensure it. Finally, the reasoning of § !2.1.2l can be made again with the more complicated non-linear 
equation, so that linearity is obtained in this case also 10 . 

This result is another illustration of the general property mentioned in the preceding section: 
here, it is the self-consistent character of the calculation that avoids the IR divergences - the 
solution automatically adjusts its low k behavior - and linearity is indeed obtained. 

2.2.3 Calculation of the linear coefficient c 

Even when linearity is proved, a difficult problem remains: the calculation of the linear coefficient 
c. Analytically, this is an intricate problem, because of its highly non-perturbative character. 
In the context of non-linear self-consistent equations, we have seen that the problem is non- 
perturbative at two levels: the choice of the non-linear equation, and the resolution of the equation 
by successive approximations. In usual perturbation theory, infrared divergences occur in all the 
integrals contained in the expression of Ek, preventing any expansion in the small parameter a/A 
around the ideal gas spectrum iy(k) = ek- When one includes higher and higher order in a/A 
terms in the energy, these IR divergences become more and more severe, so that higher powers 
of a/A in the denominator compensate exactly those in the numerator. As a result, all terms are 
comparable - see the more detailed discussion of refs. 011], which show that diagrams of all orders 
in a/A may contribute comparable amounts to the small momentum part of the spectrum, and 
therefore to the critical temperature shift. The complexity of the general problem is illustrated by 
the second leading correction: instead of being proportional to a 2 , as one could naively expect, it 
is proportional to a 2 log a, and therefore manifestly non-analytic |51 120| . 

However, if we are only interested in the leading order (linear) shift of T c , a simplification is 

that factors [ e ' 3 ( ek ~' 1 + Efc ) — l] can systematically be replaced by [/3(ek — ^ + Efc)] -1 ; in other 
words, the problem reduces to classical field theory, a case in which the temperature shift is exactly 
proportional to a to all orders 0]. This simplification can be regarded as the generalization to all 
higher order contributions of the approximation of 123() by Q24p. In the absence of non- linear a 
terms, numerical calculations for finite values of a are in principle easier, since a sometimes delicate 
extrapolation to zero a values becomes unnecessary. Nevertheless, in practice, space discretization 

10 In our proof, we have taken for granted the existence of a self-consistent expression for the effective energies, 
as in Green's function theory. The other assumption is the existence of a solution to the non-linear equation. 
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introduces another length, the lattice parameter, which may re-introduce non-linearities; a careful 
extrapolation of the lattice parameter to zero is then necessary. Taking this into account, numerical 
lattice calculations ^] E| have indeed provided a precise value, c ~ 1.3, which appears to be 
the best determination to date of this coefficient. On the other hand, classical field theory does 
not offer an analytic solution, since it is a priori not possible to select a class of well defined 
perturbation diagrams that are sufficient to give a reasonable value for c. 

Finally, another approach to the problem is given by the "large N" expansion, where one 
studies the critical temperature of a gas containing particles with many internal states, resulting 
in an order parameter with iV components; the problem is exactly soluble in the limit N — > oo 
|2"T] and the calculation of various 1/N P corrections [2j][22 allows one to extrapolate to the case 
where the particles have only one internal state (N = 2). Other perturbative- variational methods 
have also been used, see 131 1251 1271125] . 

2.3 Discussion 

We have already mentioned in § 12.11 that a spatial rearrangement of atoms with low energies 
takes place, a similar effect to that discussed by Lamb et al. ^3], but in the space of the relative 
positions of the particles (correlations) instead of their ordinary positions. The effect is maximal at 
k = 0, so that the atoms in this level condense more easily than they would in a simple repulsive 
mean field with no k dependence; of course, atoms with very low k undergo a similar effect, which 
vanishes progressively when their momentum increases beyond k c . The resulting spectrum ek + Sk 
is shown schematically in figure 5-a, with the origin of the energies at So (the value of the chemical 
potential at transition); for comparison, a purely parabolic spectrum is also shown (dashed line). 
Figure 5-b shows the corresponding variations of the populations pk, multiplied by k 2 in order 
to give directly their contribution to the total population N. At condensation, when k — > 0, the 
function goes to a constant for the ideal gas (dashed line), but to zero for the ideal gas (full line) - 
in the presence of interactions, the delta function corresponding to the presence of a Bose-Einstein 
condensate appears with no pedestal. The comparison between the situations for the ideal gas and 
the repulsive gas immediately shows that every k ^ level is less populated at the critical point 
in the presence of interactions. Therefore, An c is negative at constant T, or conversely AT C is 
positive at constant n . 

2.3.1 Role of the spectrum 

As in 012], we introduce the function U(k) by: 

U(k) = |J [S k - S ] = (2t:X 2 )- 1 [W(k) - /?e k ] (31) 

If we replace in (|10|) pk by the long wavelength approximate form [fi (ek + Sk — So)] \ we can 
express the change of the critical density created by the interactions as: 

A "^-i/*ra> (32) 

A priori, one could imagine that the value of An c depends crucially on the details of the spectrum 
near k = 0. Actually, the density change depends mostly on one parameter, the crossover point 
where U(k) becomes equal to k 2 (the point at which the energy shift introduced by the interaction 
equals the kinetic energy); it is relatively insensitive to the details of the variations of the spectrum. 
For very low values of k (<C k c ), the role of the function U(k) is to harden the free particle spectrum 
so that integrals such as (|15fl converge; therefore U(k) ^> k 2 and the integrand in l|32|l is almost 
unity. For a value k c of k comparable to a/ X 2 , see equation (I20|) , there is a crossover between the 
two functions, and one can expect that one rapidly reaches the opposite regime where k 2 3> U(k) 

1 1 Mathematically, a negative value is not totally excluded, but it would require complicated compensation effects 
and somewhat unplausible k variations of S^. 
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Figure 5: Figure (a): plot of the function ek + £k (full line), with the origin of the energies at 
So, and of the same function for the ideal gas (dashed line). Figure (b) shows the populations pk, 
multiplied by fc 2 , far above Bose Einstein condensation (curve (i)), and exactly at the condensation 
point (curve (ii)); again, the ideal gas curve is shown by a dashed line. Because of the modified 
spectrum, the curve goes to zero at the origin in the presence of repulsive interactions, unlike the 
ideal gas. 



and the integrand in l|32|l is almost zero. Altogether, we almost find an "all or nothing" regime 
where the population depletion jumps rapidly from complete to zero, so that: 

2 

An c ~ —k c (33) 

TTA Z 

which depends only on the value of k at the crossover and not on the details of the variations 
of U(k). This explains in particular why the calculation of the critical density, or of the criti- 
cal temperature, is relatively insensitive to the exact value of the universal coefficient 77 which 
characterizes the power dependence of U(k) ~ k 2 ~ v at low k. 

2.3.2 Role of exchange cycles 

Feynman has emphasized the role of exchange cycles in Bose-Einstein condensation jJU EOj : con- 
densation takes place when a macroscopic exchange cycle can propagate across the whole sample 
- see also Elser's work [ST], as well as |22 for a review. The notion of exchange cycles arises in first 
quantization, where identical particles are numbered, and where a symmetrization operator S is 
applied in a second step to obtain the correct physical states. This operator can be decomposed 
into cycles of particles, usually visualized as a closed path joining all the particles involved |32| . 
The probability of occurrence of a given path decreases rapidly when the distance between two 
consecutive particles much exceeds the thermal wavelength A; this provides the length scale for 
the maximum jump of a cycle between particles 

Reference discusses how the shift of the critical temperature introduced by the interactions 
in a homogeneous system can be interpreted in terms of exchange cycles. In an ideal Bose gas just 
above the transition point, the density fluctuations are large because exchange effects tend to bunch 
the particles together. There will be regions of space where the density is lower than average and 
where, because of the maximum length A of the jump of an exchange path, cycles will not easily 
propagate; clearly, this effect will oppose the appearance of Bose-Einstein condensation. On the 
other hand, in a dilute repulsive gas, mutual repulsions create a restoring force which tends to make 
the system more homogeneous; the occurrence of large regions of low density fluctuations becomes 
unlikely This will facilitate the propagation of large exchange cycles across the whole system, i.e., 
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Bose-Einstein condensation. Therefore, a lower density (or conversely a higher temperature) are 
needed to reach condensation 12 . 

At this point, one may wonder what is the relation between this physical explanation and 
the calculations presented above: what is relevant is the effects of IR divergences, and therefore 
physical properties of the system at very large distances; in this section, what matters is the 
correlations of particles at a distance of the order of A. This apparent paradox is removed if 
one remembers that, in our preceding discussion, we were studying one-body properties; here we 
are dealing with the propagation of exchange cycles, which depends on two-body properties. A 
phenomenon which appears at long distances for single particles can translate into shorter range 
properties for correlations. Here, we will not give a precise discussion of this point, but just develop 
a plausibility argument. We use a simplified model where the two-body density operator pjj in 
the gas is given by: 



PII {\, 2) = Ml) <g> PI (2)} [1 + P ox ] + .... (34) 

where pi is the one-body density operator and P cx is the exchange operator between particles 
1 and 2. This relation is exact for the ideal gas, but only an approximation for the interacting 
gas. Even if some effects of the interactions are already contained in (f 34|) through changes of the 
Pi's, it is clear that more terms should be added to the right of (|34l) to get an exact formula - in 
particular terms that introduce short range correlations between the particles. This approximation 
is sufficient for our present qualitative discussion. 

The direct term in the right side of (|34(l is simply a product containing no spatial correlation, 
so that we will ignore it and focus on the exchange term containing P CK with diagonal elements: 

< n,ra | Pfi(l,2) | n,r 2 > = <n,r 2 | Pl (l) ® P i{2) | r 2 ,n > (35) 

Since one-particle density operators are diagonal in the momentum representation (translational 
invariance), this term is proportional to: 

fd 3 h J d 3 fc 2 p kl Pk 2 e l ( kl - k2 Hri-r 2 ) (36) 

Therefore the two-body correlation function is independent of the sum ri + r 2 (position of the 
center of mass of the two particles), while its ri — r 2 = r dependence is given by the Fourier 
transform of the integral: 

J d 3 K P(K+k)/2 X P(K-k)/2 (37) 

which, using parity, can be expressed as a convolution integral: 



/ 



d 3 K P( k +K)/2 X /?( k -K)/2 (38) 



At the critical point, we have seen above that the distribution function p^ is the sum of 
two components: a Bose-Einstein distribution obtained within a mean-field approximation with 
a slightly negative value of the effective chemical potential (proportional to the square of the 
scattering length a), and a critical perturbation introduced for small k (comparable to k c ) by the 
correlation effects; the latter is difficult to calculate, but small. We therefore have: 

Pk = p^ + (39) 

where the first term has a width comparable to 1/A while the second has a width k c . The two-body 
correlation function now appears as the sum of three terms: a term in 



(o) 2 



a crossed term in 



12 This is actually an interesting situation: in most cases, repulsive interactions between particles tend to mask 
quantum effects, especially exchange effects (for instance, hard core potential effects reduce exchange effects in 
liquid helium three and four); for the dilute system here, they actually enhance them. 
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p k x Apt on which we will focus our interest, and finally a term in [Ap k ] , which we will ignore 
since it is smaller. 

The first term is well-known it corresponds essentially to an ideal Bose gas close to the 
transition point and contains the usual "exchange bump", with a width comparable to A, as well 
as a long exponential tail 13 . The second term is the crossed term; its Fourier transform appears 
as the convolution of the function p^j 2 by Ap k / 2 - Since Ap k / 2 is much narrower than /0 k / 2 ' the 
function Ap k / 2 may be approximated by a delta function at the origin, with a positive weight d. 
We then obtain for the Fourier transform of this term: 

d x p k/2 (40) 

This corresponds to a change of the spatial correlation function that is positive; it has basically 
the same spatial dependence as the zero order term, except that no squared function appears 
here. It therefore contains an exchange bump that extends to distances comparable to the thermal 
wavelength, but slightly further because it no longer includes a convolution of two p's; there still 
is a long exponential tail, with twice the range of the zero order tail. 

This simple model illustrates how the physics behind the reduction of the critical density may 
appear at a different scale, depending whether it is expressed in terms of single particle properties 
or in terms of correlations. In the former case, the critical phenomenon is contained in Ap k , which 
contains mostly small k's and is dominated by long distances; in the latter, Ap k disappears from 
the leading term, so that distances of the order of A remain relevant . 



3 Mean field and correlation effects combined in a trap 

As already mentioned, in the thermodynamic limit, Bosc-Einstein condensation is reached in a 
trap when the density at the center of the trap is exactly the critical density for a uniform gas, 
including of course the corrections introduced by the correlations. Two effects (the compressibility 
effect of § Q] and the correlation effect of § 0) then add to each other 14 to determine the change of 
the critical value of TV at a given temperature (or, conversely, the critical temperature at fixed N). 
Nevertheless, as remarked by Arnold and Tomasik |35|. their combination is not trivial because 
spatial effects in a trap mix up various orders in a. In this section we study the effects of the 
interactions on the density profile of a gas in a trap, assuming that the temperature is above, 
or just at the critical temperature. We start from the number density at the center, from which 
we can calculate the density profile by using an equation of state obtained within mean field 
approximation. The critical density profile is then derived by setting the density at the center to 
its critical value, including the correlation effects discussed in the preceding sections, which are 
beyond mean field theory. There is nevertheless no inconsistency in this approach; the reason is 
that the calculation of the equation of state of the gas reduces to an ordinary virial correction to 
the pressure, which is perturbative, while the calculation of the critical density remains essentially 
non perturbative. 

We will see that the specific properties of the curves of figure 1 play a role in this problem, 
in particular the large compressibility of the gas at the critical point. Roughly speaking, the 
compressibility of the gas at the center of the trap is large (~ , but remains much smaller 
in most other regions of the trap (practically independent of a) . One can then anticipate that a 
change of the pressure that is only second order in a will correspond to a change of the density 
of the same order in those regions, but first order at the center. This will allow one to meet 
the change of the condensation condition introduced by the correlations effects. But, since N is 
primarily determined by the density of the gas in those low compressibility regions, this change 
corresponds to a second order variation 15 of N. We now discuss this question more precisely. 

13 The tail can easily be obtained by using the approximation ~ f3~ 1 (e^— p) -1 , equivalent to a zero Matsubara 
frequency approximation - see for instance exercise 12.9 of 0. Since the effective chemical potential is proportional 
to a 2 at the transition point, the range of this exponential tail is of order A 2 /a. 

14 More precisely, subtract from each other, since they have opposite signs. 

15 Conversely, one can assume that N is fixed and that the temperature is reduced from above the transition; the 
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3.1 Density profile 

In a first step, we consider the number density of the gas at the center of the trap n(0) as a free 
parameter. We come back to equation Q and distinguish two cases: 
(i) Classical gas (n\ 3 -C 1) 

The relation between the chemical potential and the density is then /3/j, c g ~ ln(nA 3 ) , which 
implies J ~ 1/nA 3 , so that Q becomes: 

-V [PV + 4aA 2 n] = V ln(n) (41) 

An r integration then gives: 

n ~ exp [-(3V - 4a\ 2 n] (42) 

which is nothing but the usual Boltzmann exponential. Around the center of the trap, the density 
varies proportionally to the potential energy V: 



n{r) ~ n(0) 



1 _ PV(v) 



1 + 4aA 2 n(0) 



(if r -> 0) (43) 



(ii) Quantum non condensed gas (1 < nA 3 < C3/2) 

For simplicity, we assume that the gas is strongly degenerate at the center of the trap (n(O)A 3 
is close to C3/2)) an d we limit our study to the region of the trap where this degeneracy remains 
strong. We can then use the Mellin formula limited to lowest order (a more precise calculation 
is given in the appendix): g 3/ / 2 (exp/3/x) ~ £ 3 / 2 — 2^/—irj3^ [3IHZ!> where £3/2 is the value at the 
origin of the 33/2 function (C3/2 = ff3/2( z = 1) — 2.61..). This provides: 

J^-i-(nA 3 -n°A 3 ) 2 + .. (44) 



where: 







C3/2 A" 3 (45) 

is the critical density of the ideal gas. This result, inserted into J7J), allows us to integrate the 
gradients on each side of the equation and yields (assuming that the potential at the center of the 
trap vanishes): 



A 



This equation is second degree in n; choosing the solution which tends to n(0) when r — > 
provides: 

87ra 



n(r) = n° c + A4 



8™! 2 
n° c - n(0) + — 



%PV{v) (47) 



Equation 147|l gives the density profile of the gas as a function of r, provided that the gas is in 
the quantum regime. 

For an ideal gas (a — 0), figure 6-a shows the density profiles obtained as a function of the 
central density n(0), assuming a quadratic potential V[r) ~ r 2 ; as usual, the thermal range Rx 
is defined 16 by (r/i?^) 2 = (3V(r). If n(0) = n°, the density variations are linear in r in all the 
quantum region; for lower values of n(0), there is a parabolic variation near the center of the trap, 
in a domain that becomes larger and larger when n° — n(0) increases. 

For an interacting gas (a > 0), one can distinguish in l|47(l two possibilities, depending whether 
the potential energy j3V is larger or smaller than [(n° — n(0)) A 3 + 87ra/A] . Near the center of 

size of the atomic cloud then decreases progressively. When the system reaches a temperature close to condensation, 
the compressibility at the center becomes ~ 1/a and therefore very large, so that particles tend to accumulate more 
around this point of space than elsewhere in the trap. Consequently, a second order in a change of the temperature 
is sufficient to create a first order change of the local density. 

16 Rt is the size of the ideal classical gas in the trap at temperature T. 



16 




Figure 6: Figure (a) shows the modification of the density profile for the trapped ideal gas when 
increasing the central density until it reaches the condensation condition (full line: n(0)A 3 = £3/2! 
dashed line: n(0)A 3 = 2C 3 / 2 /3; dash-double dot line: n(0)A 3 = C3/2/2; dash-dot line: n(0)A 3 = 
£3/2/3). Figure (b) shows the corresponding density profile for the trapped repulsive gas calculated 
within mean field theory. The asymptotes at the center of the trap are calculated with the lowest 
order approximation of the Mcllin formula, equation l|47|) . 



the trap, V is small and the density variation is proportional to V (quadratic in r for a harmonic 
potential): 

- "(0) - ,3.o .0 mm + - ( if r "> °) (48) 



A 3 [0° - 7i(0))A 3 + (87ra/A)] 



This formula is similar to l)43|) but predicts significantly different results. For instance, if a is 
sufficiently small and if the central density is sufficiently close to its critical value, the change of 
the density induced by the potential becomes arbitrarily large - while it remains constant for a 
classical gas; this is a consequence of the large compressibility of the quantum gas at the origin. 
Further from the center of the trap, V is large, and (|47J) becomes: 



n(r) ~ n° + — - 



47r/?y(r) 



A (i 



[(n°-n(0))A 3 + (87ra/A)]' 



(49) 



which now predicts a square root dependence in V. 
behaviors of the density. 



Figure 6-b illustrates these two different 



3.2 Critical density profile 

We are now in position to discuss the effects of the interactions on the density profile at the 
transition point, at least in the central region of the trap where the gas remains strongly degenerate. 
We set the central density to its critical value n c and use the above equations to calculate n(r). 
This calculation is of course not exact, for different reasons: first, (|47|l is only valid within mean 
field theory; second, the use of the first correction only in the Mellin formula limits us to a 
sufficiently degenerate gas. Both approximations could be improved: the first by including second 
order perturbative corrections to the equation of state, as done in ref. |35j : the second by including 
more terms in the Mellin formula, as done in the appendix. Our calculation nevertheless remains 
sufficient for a qualitative discussion. 

Let us first ignore the effects of the correlations on the critical density at the center. We then 
set n(0) = n° and see that the density corrections are positive and first order in a everywhere in 
the trap - second term in the right side of eq. I|49(l - except in a small region around the center 
where the density change is still positive but oc 1/a - fraction in the right side of eq. (|48[l. For 
a harmonic potential, this small region has a size proportional to a, so that the corresponding 
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contribution to N is only second order in a. Therefore, the main contribution arises from all the 
rest of the trap and provides a first order in a change of N (at constant temperature) ; conversely, 
calculating the temperature shift at constant N provides JHJ, as shown in ref. 

We now take into account the effects of the correlations on the critical density at the center 
of the trap. We then need to choose a slightly smaller value of n(0) since, according to (|26|l . n c 
is reduced by an amount that is first order in a. With respect to the calculation of the preceding 
paragraph, eq. Q48JI predicts an additional negative change, first order in a, of the density around 
the center of the trap. On the other hand, eq. (|49|) predicts only a second order change induced 
by the new value of n(0); the crossover between the regions of space where these equations apply 
is obtained when the potential energy is second order in a, i.e., at a first order distance from the 
origin if the potential is harmonic. This corresponds only to a fourth order correction to the total 
number of particles N, because the volume of this region is proportional to a 3 for a harmonic 
potential; this correction remains therefore negligible when compared to the second order in a 
contribution from all other regions in the trap. Altogether, we see that the additional change of 
N introduced by the correlation effects is therefore only second order in a; this is to be compared 
with the first order correction introduced by mean field effects contained in the term in a/VA 4 on 
the right side of 

We therefore recover the conclusions of Arnold and Tomasik [123 : the two N changes subtract 
from each other, but not at the same order in a. This offset of one order illustrates again the 
relatively loose connection between the total number of atoms in the trap N and the physics at 
the center, where Bose-Einstein first takes place: most of the particles contributing to N are too 
far from condensation to play a role in it. If the central density n(0) were accessible, it would 
provide a more adequate variable than N for the study of the BEC phenomenon. See ref. for 
a proposal of a method (adiabatic ramping down of the trap frequency) aimed at enhancing the 
visibility of the correlation effects against the background of mean field effects. 

4 Conclusion 

The calculation of the critical temperature of a dilute Bose gas has several distinguishing features. 
One could of course try to go further than the leading a correction to the critical temperature and 
attempt to describe analytically the full density variations obtained by path integral Monte Carlo 
methods in [33]. Quantitatively it is clear that at high densities the blocking of exchange created 
by the repulsive hard core part of the potential tends to reduce the critical temperature, an effect 
which goes in the opposite direction to what happens for dilute gases. Ref. {38) gives a discussion 
of this effect in terms of change of the effective mass of the particles, arising from a k dependence 
of the exchange term in the expression of the mean field. It would be interesting to make this idea 
more quantitative and to explore the specific details of BEC in all density regimes. 
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APPENDIX 

A better approximation for g 3 / 2 than that used in § [3] is: 

53 /2 (exp/?/i) ~ C3/2 - 2y/-n/3fi - b f3(i + ... (50) 

with b = 1.46.. [2 EH- Adding a linear term to the square root actually provides values which 
are reasonably accurate until f3/j, reaches values where the classical regime is obtained < — 1). 
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Figure 7: Full line: mean field density profile calculated exactly (equation (@J). Dashed line: same 
profile with the lowest order approximation of the Mellin formula (equation (|47(0 . Dotted line: 
same profile with the inclusion of the linear correction (equation 



Inverting the relation nA 3 = g 3 / 2 (exp /?/j, G ff) then gives: 



-PfJ-eS - 



2 + h - (n - n° c ) A 3 - 2^1 + | (n - A3 



(51) 



which defines the function J(nA 3 ). We therefore have: 



J 



^1 + 6A3 (n-n° c )/TT 
Relation J7J then becomes, after integration of the gradients: 
-0V = 4aA 2 (n - n(0)) 



(52) 



(n - n(0)) - + 6A3 (n - ^) /tt + ^^1 + 6A3 (n(0) - n°) /tt 



(53) 



This equation provides a direct relation between the potential energy and the local density. For 
instance, for a quadratic potential, the position r is obtained by taking the square root of the right 
side, multiplied by some coefficient. This allows one to make a precise calculation of the density 
profile in all regions of the trap where the gas is in a quantum regime. Figure 7 gives a comparison 
between the density profile obtained in the mean field approximation with l|53|l and (|47|l . 
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